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The paper is an attempt to relate two vast areas of the applicability of the renormalization 
group (RG): field theoretic models and partial differential equations. It is shown that the Green 
function of a nonlinear diffusion equation can be viewed as a correlation function in a field-theoretic 
model with an ultralocal term, concentrated at a spacetime point. This field theory is shown to be 
multiplicatively renormalizable, so that the RG equations can be derived in a standard fashion, and 
the RG functions (the j3 function and anomalous dimensions) can be calculated within a controlled 

■ approximation. A direct calculation carried out in the two-loop approximation for the nonlinearity 
of the form (f> a , where a > 1 is not necessarily integer, confirms the validity and self-consistency of 
the approach. The explicit self-similar solution is obtained for the infrared asymptotic region, with 
exactly known exponents; its range of validity and relationship to previous treatments are briefly 
discussed. 

(N ■ PACS numbers: 05.10.Cc, ll.IO.Gh, 02.30. Jr 
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■ I. INTRODUCTION 

d . 

The renormalization group (RG) has proved to be the most efficient tool for studying self-similar scaling behavior. 
First appeared within the context of quantum field theory Jjj] , it was then successfully applied to a variety of problems 
as disparate as phase transitions, polymer dilutes, random walks, hydrodynamical turbulence, growth processes, and 
so on; see, e.g., the monorgaphs ^|J^], the proceedings Q, and references therein. 

The most powerful and well-developed formulation of the RG is the field theoretic one; see (|-[§. It is this version of 
\Q , the RG that is simplest and most convenient in practical calculations, especially in higher orders. It is also important 

■ that it has a reliable basis in the form of quantum-field renormalization theory, including the renormalization of 
composite operators and operator product expansion. For this reason, the first step in the RG analysis of a given 
problem is to reformulate it as a field theoretic model. This means that the quantities under study should be 
represented as functional averages with the weight exp 5 (</>), where <f> is a classical random field (or set of fields) 
and S((f>) is certain action functional. For parabolic differential equations with an additive random source, such a 
formulation is provided by the well-known Martin-Siggia-Rose formalism, see J5]fl . In problems involving fluctuation 
effects in chemical reactions the somewhat more complicated approach of Doi |jj (see also |||]]) nas a l so been widely 
used j[o[[l2). No general recipe, however, seems to exist to cast a nonlinear problem to a field-theoretic form. 

Such a reformulation, however, is by no means superfluous: once the field theoretic formulation has been found, it 
becomes possible to apply standard tools (power counting of the 1-irreducible correlation functions etc) to verify the 
renormalizability of the model, i.e., the applicability of the RG technique, to derive corresponding RG equations, and to 
calculate its coefficients (beta functions and anomalous dimensions) within controlled approximations. An instructive 
example is provided by the model of the so-called true self-avoiding random walks [p^|-p^[. After its field theoretic 
formulation had been found [fblf , it became clear that the model in its original formulation was not renormalizable, 
and the direct application of the RG to it would lead to completely erroneous results. The renormalizable version of 
the model can be obtained by adding of infinitely many terms to the original action; see Jig] . 

It has long been known, however, that symmetries of the RG type also appear in various physical problems described 
by ordinary or partial differential equations and intcgro-diffcrcntial equations, whose solutions exhibit self-similar 
scaling behavior [|l6|. Since then, the list of such problems has been essentially increased; see, e.g., an< ^ 
references therein. As a rule, the field theoretic formulation for these models does not exist (or, at least, is not 
known), and derivation of the corresponding RG equations is a nontrivial task. Quoting the authors of |24[ , "the 
procedure of revealing RG transformations ... in any partial case ... up to now is not a regular one. In practice, 
it needs some imagination and atypical manipulation 'invented' for every particular case." In Ref. a general 
approach was proposed to constructing RG symmetries for certain classes of partial differential equations, but its 
relationship to the field theoretic RG techniques is not clear. 

The present paper is an attempt to 'bridge the gap' between these two vast areas of the applicability of the RG: 
field theoretic models and partial differential equations. To be specific, we shall consider nonlinear diffusion equation 
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of the form 



d t ct> = v d 2 (t> + v{(t>), 



(i.i) 



where 4>{x) = </>(£, x) is a scalar field, 2/0 is the diffusion coefficient, <9 2 is the Laplace operator, and V(4>) is some 
nonlinearity dependent on the field <j) and its spatial derivatives. Within the RG context, various special examples of 
Eq. (1.1) were studied earlier in [jl8] p2| . In practical calculations, we shall confine ourselves to the nonlinearity of 
the form V(<fi) — — Ao0 Q , where a > 1 is not necessarily integer. 

We shall show that the problem (1.1) can be cast into a field theoretic model and apply the standard RG formalism 
to it to establish the scaling behavior and to calculate corresponding anomalous dimensions. Then we shall discuss 
the range of applicability of the results obtained and their relationship to the previous RG treatments of the model. 



II. FIELD THEORETIC FORMULATION AND RENORMALIZATION OF THE PROBLEM 

We begin the analysis of the Cauchy problem ( p_.l| ) with a localized initial condition which corresponds to the 
equation 

d t G = u d 2 G + V(G) + S(x - x ) (2.1) 
for the Green function G(a;|eo). It will be shown later that the large-scale asymptotic behavior of this problem 



survives for all integrable initial conditions (i.e. such that J<ix0(O,x) converges). In Eq. (2.1) we denote 8{x — xq) 



5(t — to)S^(x — xq), where d is the dimensionality of the x space, and gq = {xq, vq, Aq} is the full set of pa ram eters 



The functional derivation of the MSR formalism |6| can be adopted to represent the solution of Eq. (2J) as a 
functional integral over the doubled set of fields, <fi and 4>': 



G(x\e a ) = J V4>' I Vcp <j>(x) exp[S(<f/, <t>) + 4>'(x Q )] . (2.2) 
Here the normalization constant is included into the differential 2?</>'2?0, the action functional has the form 

S (</>', 0) = fdx cf>'{x)[-d t ci>{x) + v Q d 2 4(x) + V (4>(x)) }, (2.3) 



with dx = dtdx. The last term in (2.1) can be treated as an addition to the 'interaction' V{4>) and gives rise to the 



last term in the exponential of Eq. (|2.2|). The t erm quadratic in <//, typical to the MSR actions, is absent in (|2.3|) 



owing to the absence of the random force in Eq. (1.1) 



Representation (|2.2|) shows that the Gree n fu nction (2.1) can be viewed as the correlation function (4>(x) exp (j)'(xo)) 



in the field-theoretic model with the action (2.3). It is not convenient, however, to deal with the exponential composite 
operator exp</>'. A more useful interpretation is the following: the integral ( |2.2| ) describes the correlation function 
(4>(x)} for the extended action S' = S + (f>'(xo) with an 'ultralocaP interaction term concentrated on a single spacetime 
point Xq. 

The renormalization of field theoretic models with ultralocal terms, concentrated on surfaces, was studied in Ref. 
p5[ in detail. The analysis of |2^], which we also naturally take to apply to our case, has shown that the standard 
renormalization theory is applicable to such models, with some obvious modification (see below). 

The analysis of ult ravi olet (UV) divergences is based on the analysis of canonical dimensions; see Dynamical 
models of the type ( |2.3| ), in contrast to static models, have two scales, the length scale L and the time scale T. 
Therefore, the canonical dimension of any quantity F (a field or a parameter in the action functional) is described by 
two numbers, the momentum dimension dp and the frequency dimension d%, determined so that [F] - [L\~ d ' F [T]~ d ^\ 
The dimensions are found from the obvious normalization conditions d\ = — <i£ = 1, d% = d^ = 0, d£, = d^ = 0, 
d^ — — d^ = 1, and from the requirement that each term of the action functional be dimensionless (with respect to the 
momentum and frequency dimensions separately). Then, based on dp and d F , one can introduce the total canonical 
dimension dF = d k F + 2d F (in the free theory, dt oc d 2 ), which plays in the theory of renormalization of dynamical 
models the same role as the conventional (momentum) dimension does in static problems; see ^,|). 

Now let us turn to the special case of the model (f2.1|) with the extended action of the form 



S'{4>\4>)= f dx<t>'(x){-d t cl>(x)+i,od 2 cP(x)-gov <t> a (x)}+<p'(x ), (2.4) 
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where we have introduced the new parameter go = Aq/Vo, which plays the part of the coup ling constant (a formal 



small parameter of the ordinary perturbation theory). Canonical dimensions for the model (2^) are given in Table 
H including the dimensions of renormalized parameters, which will appear later on. From Table [| it follows that 
the model is logarithmic (th e coupling constant go is dimensionless) for 2 + d(l — a) = 0. In what follows, we fix 



the exponent a in Eq. (2.4) and consider the model in variable space dimension d = (2 — s)/(a — 1). Then the 
UV divergences take on the form of the poles in e = 2 + d(l — a) in the correlation functions. The 'interaction' is 
therefore irrelevant (in the sense of Wilson) for e < 0, marginal (logarithmic) for e — 0, and relevant for e > 0; cf. 
the analysis in Ref. [fl9|| . This means that for s > 0, the ordinary perturbation expansion (i.e., series in go) fails to 
give correct infrared (IR) behavior and has to be summed up. The desired summation can be accomplished using the 
renormalization group. 

It is common wisdom of the renormalization theory that for the analysis of UV divergences of all correlation 
functions of the fields <fi and cf>' it is sufficient to consider one-particle-irreducible (1PI) correlation functions, whose 
graphical representation contains only graphs which remain connected after removal of one (arbitrary) line (i.e. a 
free-field correlation or response function) of the graph. 

The total canonical dimension of an arbitrary 1PI correlation function 

T( Xl , ■ ■ • , XN , y u ■ ■ ■ , y N „ Xo ) = 5<p{xi) ... H{xN)HI[yi) ... HI[yNi) . (2-5) 

in the time-coordinate representation is given by the relation 

d r = N(d + 2-d 4> )+N'(d + 2~d 4> ,), (2.6) 

where N and N' are the numbers of corresponding fields. In (^1]) T((/>, tfi') is the (dimensionless) generating functional 
of 1PI Green functions. It should be noted, however, that due to the presence of the ultralocal term in the action, 
the functional T((j), cf)') is not the Legendre transform of the functional W(J,J') = lnQ(J, J'), where Q(J : J') — 
J "Dtp' J T>cf> exp [S" ((/)', <j>) + Jcf> + J' </>')] is the generating functional of Green functions of the model. Moreover, 
contrary to the usual field theories, the functional lnQ(J, J') does not include all connected graphs of G{J, J')- By 
definition of the generating functional, the 1PI Green function with N external cf> legs and N' external (/>' legs may be 
obtained by N functional differentiations of T(<p, <p') with respect to the field <fi and N' differentiations with respect to 
4>' . The canonical dimensions of the functional derivatives are related to the dimensions of the corresponding fields as 
d k [8/6(j)] = d — dh, d u [5/8(1)] = 1 — tf^, and similarly for the auxiliary field <f>' . Then the total canonical dimension of 
the function ( |2.5| ) in the frequency-momentum representation (obtained by the Fourier transformation with respect 



to all N + N' independent differences of the time and coordinate arguments) is obtained from (2.6) by subtracting 
the term (N + N')(d + 2) and has the form 

d r = -d^N - d^N' = -dN, (2.7) 

where the data fr om Table § are used in the last equality. 



The quantity (2.7) is the formal index of the UV divergence for the function T. Like for usual (local) models, 
superficial UV divergences, whose removal requires counterterms, can be present only in those functions T for which 
S = dr\ e= o is a non-negative integer; see 

From Eq. (|]t]) we conclude that for any positive d, such divergences can exist only in the 1PI functions with N = 
and arbitrary value of N'. For all these functions 6 = 0, that is, the divergences are logarithmic and the corresponding 
counterterms in the frequency-momentum representation are constants. 



At first glance, we have established that the model (2.4) requires infinitely many counterterms, and hence it is not 



renormalizable. However, it turns out to be sufficient to renormalize the 1PI Green function T(x;xo) only to render 
the model finite, as we shall now show. 

The first few Feynman diagrams of G are shown in Fig. I for a = 2; the symmetry coefficients are shown for general 
a (it would be embarrassing to depict the diagrams for fractional a, but the idea is the same). The lines with a slash 
denote the bare propagator 

AU r.\ IAAJ\ exp[-r 2 /4^ t] 

A(M) = (#)o= (w)d/2 , (2.8) 

the end with a slash corresponds to the field <f>' , and the end without a slash corresponds to <f>. The initial (left) point 
in each diagram corresponds to x, and the final (right) point with variable number of attached lines corresponds to 
Xq. The crucial point is that, as is easily seen from Fig. I, all possible 1PI sub-diagrams entering into the diagrams 
of G belong to the only 1PI function T{x;xo); no other 1PI functions are involved. The function G appears to be 
'closed with respect to the renormalization,' i.e., we can eliminate their UV divergences by the only counterterm 
corresponding to its 1PI part r(a;;xo). 
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FIG. 1. First Feynman graphs for the Green function of the nonlinear diffusion equation (0) for V(G) = -AoG 2 . 

Moreover, the renormalization of the only function T(x; Xq) is in fact sufficient to completely renormalize all functions 
with N' > 1. A typical diagram for N' — 3 is shown in Fig. II. It is clear that any such diagram reduces to a product 
of blocks that belong to the simplest function with N' = 1 (we recall that there is no integration over xo, the only 
point that connects the blocks). Therefore, the diagram contains no superficial divergences; all its divergences are 
those of the subdiagrams and they are completely removed by the renormalization of the function with N' = 1. This 
is equally true for any diagram of any function with N' > 1. 



FIG. 2. A three-loop Feynman graph for the three-point correlation function of the nonlinear diffusion equation (2.1) for 
V(G) — —AoG 2 illustrating the factorization property T(xi, X2, £3; xo) — T(xi;xo)T(x2;xo)r(x3;xo). 

In the generic case all the loops are created by the presence of a single local vertex with any number of <j)' legs, from 
which continous chains of retarded diffusion propagators emanate. Due to the structure of the nonlinear term these 
chains do not branch, but they may merge (the single (//-field in the nonlinearity allows only one outgoing propagator 
from each ordinary vertex, whereas up to a incoming chains are allowed) . A little reflection along the lines sketched 
above shows then that all divergent 1PI Green functions are factorized: 

T(xx, ■ ■ ■ ,x N ;x ) = r(xi,x ) ■ ■■T(x N ;x ) . 

Thus, we are left with the only counterterm to the function T(x]Xq). It is constant (see above), which in the time- 
coordinate representation corresponds to the function 5(x — x ) = 5(t — t )5 t - d ' > (x — x ). In the action functional, 
after the integration over the field argument, this gives 4>'{xq). Such term is present in the extended action ( |2.4| ). 
so that our model is renormalized multiplicatively, with the only renormalization constant, which we denote Z . The 
renormalized action has the form 

S' R [<f/,<f>)= J dx cb'{x)^-d t (b(x) + ^d 2 (b(x)-giy f i e cl) a (x)'j+Z(j) , (xo). (2.9) 

Here and below the g and v are the renormalized analogs of the bare parameters, \x is the reference mass in the minimal 
subtraction (MS) scheme, which we use in practical calculations, and the constant Z depends on the dimensionlcss 
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para meters g, a and e. The renormalized Green function G R , which is finite for e — > 0, is given by the representation 
(U) with the substitution 5" -> S' R . 

If we now replace the local initial condition with an integrable one: 0(0, x) = a(x), then - after Fourier transforming 
- we obtain wave-vector integrals in which all the propagator lines starting from the initial condition contain a 
multiplicative factor 5(k) = J 5x e - a(x). For the large-scale asymptotic analysis using RG it is sufficient to keep 
the leading small wave-number terms in all the lines which amounts to the replacement 5(k) — > 5(0) = Jrfxa(x), and 
we thus return to loop integrals of the problem with localized initial condition in which Jrfxa(x) is the amplitude of 
the initial S function. 

To clarify the idea, consider the one-loop graph of Fig. I, whose analytic expression with the initial condition 
0(0, x) = a(x) is 

oo 

T« (t, x) = -A fdyj dt'A(t - t', x - y) J d yi A(t', y - yi)a(yi) J dy 2 A(t', y - y 2 )a(y 2 ) . 
o 

Here, A(i,x) is the diffusion propagator ( |2.8| ). Fourier transforming T^(t,x) with respect to x we arrive at the 
expression 



oo 

r«(t,k) = -A jdt'A(t-t',k) y^ 3 A(t / 1 q)5(q)A(f > k-q)a(k-q) ! 



(2.10) 



where A(i,k) is the spatial Fourier transform of the diffusion kernel (2J3). From the point of view of RG, the IR 



relevant terms are give n by the leading terms of the gradient expansion of the initial condition: 5(q) = 5(0) + o(q)/q. 



This allows to replace ( 2.10 ) by 

oo 

r«(t,k)~-Ao Jdt'A(t-t',k) |^_A(t',q)A(t',k-q)5 2 (0) 



which corresponds to the localized initial condition with the amplitude 5(0) = J"5xa(x). 



III. RG EQUATIONS AND RG FUNCTIONS 



S'(Z<j>',Z- 1 cf ) ,eo) = S' R 



It follows from Eqs. (|2.3| ), ([2.4|), and (2J3) that the original and renormalized action functionals satisfy the relation 



,e,fj,), if the bare and renormalized parameters are related as follows: 



9a = 9^ e Z 



€ rya — 1 



(3.1) 



with the only renormalization constant Z from Eq. (pj|). This implies the relation G(eo) = Z 1 Gfl(e,/i) for the 



corresponding Green functions in Eq. (2.2); i.e., this quantity is multiplicatively renormalizable. We use T>^ to denote 
the differential operation fid^ for fixed eo and operate on both sides of this equation with it. This gives the basic RG 
equation: 



V^+(3(g)d g -l{g) G R (e,f,)=0, 



(3.2) 



where T>^ + P(g)d g is nothing else than the operation T>^ expressed in the renormalized variables. In Eq. (3.2), we 
have written T> x = xd x for any variable x, and the RG functions (the (3 function and the anomalous dimensions 7) 
are defined as 



7(5) 



V^lnZ, 



fig = 



(a - 1)7(5) 



(3.3) 



The relation between (3 and 7 results from the definitions and the relations ( p.l| ) . 

We shall see below that, for small e > 0, an IR stable fixed point 3* of the RG equation (3.2) exists in the physical 
region g > i.e., /?(<?*) = 0, (3'{g*) > 0. The functions G and G R coincide up to a const ant(i.e., independent of the 
time and space variables) factor Z and the choice of the parameters (bare eo or renormalized e, /x) and can equally 
be used in the analysis of the IR behavior. The general solution of the RG equations is discussed in detail, e.g., in 
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[§,!)■ It follows from this solution tha t, w hen an IR stable fixed point is present, the leading term of the IR behavior 
of the function Gr oc G satisfies Eq. (3.2) with the substitution g — > g*: 



Gfl(e,/x) = 0. 



(3.4) 



In our case, the value of the anomalous dimension at the fixed point is found exactly owing to the relation between 
P and 7 in Eq. (O): 



7 * = 7 (.9*) - -e/(a - 1) = d - 2/ (a - 1). 



(3.5) 



Dimensional considerations yield Gn(t,r) = (vt) d ^ 2 £ > (l/t/i 2 iy, r 2 /tv), where £ is some function of dimensionlcss 
variables. The dependence on g is not displayed explicitly, because the derivatives with respect to this parameter do not 

enter into Eq. (3.4). It follows from Eq. ( 3.4 ) that £ satisfies — at the fixed point — the equation T) s — r y*/2 £(s,y) = 



0, its gen eral solution is £(s, y) 



- «7*/2 



x(y)i where X is an arbitrary function of the second variable y. For the Green 



function (2J2) we then obtain 

G(t,r) ~ G R (t,r) ~ t^' 2 ^' '' 2 X {r 2 jtv) = r^^r 2 /tv), 

where the form on the 'scaling function' xi 7 " 2 l^ v ) i s n °t determined by the equation ( |3.4[ ). 
parameters u, /i can be easily restored from the dimensionality considerations (see Table I): 

G(t,r) ~ (vot)- 1 ^- 1 ) X (r 2 /tv ). 



The dependence on the 



(3.6) 



Although the value of 7* in Eq. ( pJty and the solution (3.6) have been obtained without practical calculation of the 
constant Z and functions (3.3), such calculation is needed to check the existence, positivity and IR stability of the 
fixed point. Within the e expansion, these facts can be verified already in the simplest one-loop calculation. 

In order to check the validity and self-consistency of the approach, we calculated the constant Z up to the two-loop 
approximation. The calculation is performed in the frequency-momentum (cj, k) representation and calls for the 
formulas derived in Ref. [^6| for a model of critical dynamics. 

Two key points are as follows: the convolution of two functions of the form F(a; a) = (— iuj a + k 2 )~ a is a function 
of the same form, 

F(a; a) * F(/3; b) = K(a, f3; a, b) F(a + (3- d/2 - 1; a + b) 

where a and b are both positive and the coefficient has the form: 

K(a, /3; a, b) = a d ' 2 - a b d / 2 -^{a + 6) Q+ ^- d - 1 r(a +/3-d/2- l)/r(a)r(/3)(4^) rf / 2 , 

while the product of two such functions can be represented as a single integral of a function of the same form with 
the aid of the generalized Feynman formula: 



F(a; a) ■ F(/3; b) 



ds s a - L (l - sf- 1 F{a + (3:,as + 6(1 - s)). 



r(a)r(/9) Jo 

For the sake of brevity, below we give only the final result: 



u 

e 



Z = 1 

where we have introduced a new coupling constant 



au u ~ 
2^~ Ye a 



0{u% 



_£_ a -V("-i) 



(47T) 



and have written = a In a + a a ^ a with the convergent single integral 



1 ds 



(s(c 



1+1 



(2-a)/(o-l) 



((s + l)(a 



1+1 



a -l/(c-l)\ 



(3.7) 



(3.8) 



(3.9) 



in particular, I 2 = 21n(4/3) and I 3 = 6[ln(3 - V$) + ln(3/2)] . 
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Then for the corresponding beta function we obtain (3 U 



T>,,u = —u [e + PAhiZ 01 - 1 ], where we have used the 
last relation in Eq. (3.1) and the fact that 2? M = (3 u d u for the functions dependent only on u. This yields 

-£tt 



(iu{u) = 



(3.10) 



l + (a-l)D u lnZ 
Substituting Eq. (|77j) into Eq. ( ]3.10| ) gives 

/ 3 tl (u) = -w[e-u(a-l) + u 2 (a-l)/ ct ] + 0(u 4 ). (3.11) 

Note that the poles in e in the constant Z cancel out in the function ( |3.11 ) ; this is a manifestation of the general fact 
that the RG functions must be UV finite, i.e., finite as e — > 0. The cancellation is possible by virtue of the correlation 
that exists between the u/e and (u/e) 2 terms in Eq. (3.7) and can be used as an additional check of the consistency 
of the approa ch. T he simple (linear) dependence on e is a feature specific to the MS scheme. 
From Eq. (3.11) we find an explicit expression for the coordinate of the fixed point: 



+ 0(e 3 ). 



{a -I) ' (a-1) 2 

As already said above, for small positive e and a > 1 the fixed point is positive and IR stable: (3' u (u*) 



(3.12) 
0(e 2 ). 



IV. DISCUSSION 

We have applied the field theoretic renormalization group to the non-stochastic differential equation ( |2.l| ) and 
established the scaling behavior in the IR asymptotic range, as a consequence of the existence of the IR stable fixed 
point in the physical range of parameters. The same asymptotic behavior is shown to be valid for integrable initial 
conditions which thus constitute the universality class of this fixed point. 

The key points are the formulation of the problem as a field theoretic model with an ultralocal term concentrated 
at a spacetime point and the fact that this model appears multiplicatively renormalizable, in spite of the naive power 
counting that indicates nonrenormalizability. 

The two-loop calculation confirms internal consistency of the approach. 

Simple explicit form of the scaling dimensions follows from the fact that there is only one independent renormal- 
ization constant in the problem. In particular, this explains a simple value z = 2 of the exponent in the argument 
2 jt x l z of the scaling function (3.6) (in models of dynamical critical phenomena 0,0] and some models of nonlinear 



r 



diffusion |27j this exponent differs from two). 

Recent ly, i t has been conjectured p2| that the dynamic exponent z ^ 2 in the present problem. Our asymptotic 



solution (3.6), however, does not predict any deviation from the canonical value z — 2, since there is no renormalization 
of the diffusion coefficient in the MS scheme we have used. In Ref. |2^] with the use of a different renormalization 
procedure it was concluded that z — 2 = 0(e 2 ). We think, however, that it is not consistent to prescribe physical 
quantities values of the order 0{e 2 ) on the basis of the one-loop calculation carried out in Ref. |2j|], but a two-loop 
analysis is required for this accuracy. 



The RG analysis allows one to derive the RG equation rigorously and to prove that the behavior ( p.6[ ) is indeed 
realized for e > 0, go > Q in the IR asymptotic range, specified by the relations t ~ r 2 and r <C r), where rj ~ g^ 1 ^ is 



the UV scale. The general solution of Eq. (|3.2|) interpolates between the ordinary perturbation theory for Eq. (2.1) 
and the self-similar asymptotic expression ( |3.6| ). The scaling function x(u) can be calculated within the e expansion; 
in the lowest order one easily obtains x(u) = ex P[~(?//2) 2 ] + 0(e). 

We hope that the ideas presented above might be useful in other models containing ultralocal contributions, which 
have several charges and hence richer IR behavior. Another direction of generalization would be the analysis of Green 
functions of vector quantities. 
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TABLE I. Canonical dimensions of the fields and parameters in the model 
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